Thermal state entanglement in harmonic lattices 
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We investigate the entanglement properties of thermal states of the harmonic lattice in one, two 
and three dimensions. We establish the value of the critical temperature for entanglement between 
neighbouring sites and give physical reasons. Further sites are shown to be entangled only due to 
boundary effects. Other forms of entanglement are addressed in the second part of the paper by 
using the energy as witness of entanglement. We close with a comprehensive diagram showing the 
different phases of entanglement versus complete separability and propose techniques to swap and 
tune entanglement experimentally. 



I. INTRODUCTION 

Entanglement is a purely quantum phenomenon essen- 
tial for understanding how information is shared in quan- 
tum systems and implementing new computing technolo- 
gies such as measurement-based computation. The ex- 
istence of entanglement in a composite system implies 
strong correlations between its parts; so strong that no 
classical model can explain them. When classical corre- 
lations become strong in a condensed matter system, i.e. 
when they stretch over large distances, a new 'phase' 
emerges which manifests in qualitatively different be- 
haviour of the material. Phase transitions happen when 
a parameter, such as the temperature, reaches a criti- 
cal value. In this paper we will determine the critical 
temperature below which entanglement is present and at 
which a quantum-correlated phase emerges in the har- 
monic lattice. 

The entanglement properties of the bisected one- 
dimensional harmonic chain have been analysed in [l[ 
using extensive mathematical formalism to discuss the 
ground state (see also @]). The authors also provide a 
remark on the situation for thermal states, however, the 
treatment there is rather short. In the present paper 
we illuminate the thermal case in greater detail by first 
discussing two-site entanglement and then any entangle- 
ment for finite size chains as well as in the thermody- 
namic limit, for varying external potential and in one, 
two and three dimensions. Our intention is to gain in- 
tuition about the physics behind the quantum-classical 
transition at T cr j t . We explain how the excitation and 
mixing of the phononic modes degrades the entanglement 
unless the lattice is in the intrinsically 'quantum phase' 
below T cr i t where quantum zero-point fluctuations dom- 
inate. We find that entanglement depends on external 
parameters as well as size and propose that manipulat- 
ing these will enable techniques of tuning and swapping 
entanglement in trapped ions. 

The outline of the paper is as follows. In Section HI1 the 
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harmonic lattice model is introduced and we summarise 
technical methods that build the fundament for the re- 
sults discussed later. Readers familiar with these tech- 
niques may skip to Section Hill where we derive analytical 
expressions for the entanglement of the one-dimensional 
harmonic chain in the thermodynamic limit and at zero 
temperature, and then approximate the critical temper- 
ature for nearest neighbour entanglement. In Section HV1 
we present comprehensive numerical results for two-site 
entanglement in the one-dimensional chain. The value of 
the respective critical temperatures is explained and fi- 
nite chain sizes are discussed. In Scction|V]we investigate 
the dependence of the nearest neighbour entanglement in 
the one-dimensional chain when increasing an external 
trapping potential. Finally, this behaviour is compared 
with the behaviour of entanglement in the two- and three- 
dimensional lattice. 

In Section IVII a different approach to detect entangle- 
ment is taken. Using the energy as a witness for any 
entanglement (not just between two sites) we derive the 
minimal energy bound for completely separable thermal 
states. We then establish the critical temperature for wit- 
nessed entanglement for the physical sites and find that it 
is independent of the dimension of the lattice. In Section 
IVIII we show and compare the two temperature thresh- 
olds derived in the present paper with the threshold for 
block-entanglement found in [1| and the threshold for full 
separability derived in 0] in a comprehensive entangle- 
ment phase diagram. We draw concluding remarks in 
Section IVIIII propose possible uses of our results in con- 
trolling and storing information in trapped ion systems 
and discuss the notion of an entangled phase. 



II. THE MODEL 

A. Hamiltonian and phononic spectrum 

A harmonic lattice can be realised, for example, as 
a system of ions that are trapped at N sites by a har- 
monic trapping potential 5. The ions shall couple to their 
nearest neighbours with equal strength u> throughout the 
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lattice. The Hamiltonian for the rf-dimensional lattice is 




H = 



Here is the operator for the deviation from the equi- 
librium position of each ion. The sites are indexed by 
R = na with a being the lattice constant and the com- 
ponents of n, running over n,j = 1,2,..., -\/N. Assuming 
square geometry of the lattice the length in each spatial 
direction is L = -\fN a and m is the effective mass of 
the ions, is the momentum operator corresponding 

to so that for all components i and j the commu- 

-M) Ash 



tation relations hold [u { ~',p { ij] = ihd^ ^, Si.j . The first 
sum in Eq.([T]) is taken over all N sites, and the second 
sum takes all neighbouring pairs of sites, (R, R') , that are 
harmonically coupled. To simplify the calculation we will 
adopt periodic boundary conditions and additionally re- 
quire -\fN to be an odd number for the sake of definitive 
calculations. 

The Hamiltonian can be diagonalised by applying first 
a Fourier transformation followed by a Bogoliubov trans- 
formation into a new set of position and momentum op- 
erators, see for instance The frequencies of the vi- 
brational modes of the lattice, the phonons, are given 
by 



and second moments, i.e. the expectation values of the 
canonical operators themselves and the expectation val- 
ues of any combination of two canonical operators. The 
first moments can be changed by choice of the coordinate- 
system, i.e. by applying local, single mode displacements, 
and their absolute value is of no great physical signifi- 
cance. Only their relative values obey conservation laws. 
Correlations between modes are determined by higher 
moments and for Gaussian states the second moments 
are hence vital for the discussion of entanglement. The 
second moments of a state p with modes 1, 2, .... and 
canonical operators Q = (ui,pi, U2,P2, •••) ar e collected 
in the covariance matrix, whose elements are given by 



ijfc(p) = (QjQk + QkQj) P - 2(Q j ) p (Q k ) f 



(4) 



r is a real, symmetric and positive matrix which reveals 
significant properties of the state such as the occurrence 
and amount of entanglement in the state. The advantage 
of formulating the discussion in terms of the covariance 
matrix is the substantial reduction of parameters, 2x 
(number of modes), compared to the 2( number of modos > 
for the density matrix. 

For the harmonic lattice the covariance matrix for two 
sites R and P can be found using the transformation from 
real space into phononic modes in which the Hamiltonian 
is diagonal. Then expectation values are easily calculated 
and one obtains 



2w< 



with L 



■/N-i 



■ ■ , v/ ^~ 1 for a d-dimensional square 
lattice with TV sites, each having 2d neighbours. It will 
turn out that the frequency spectrum of the phononic 
modes is key to the value of the critical temperature for 
nearest neighbour entanglement and the entanglement 
properties can be modulated by varying the interaction 
and the trapping potential. (In the following vector ar- 
rows indicating the dimension will be dropped and the 
dimension of the lattice will be stated explicitly when 
necessary.) 



B. Thermal states and covariance matrix 



A E 

B F 

E A 

F B 



(5) 



where the entries are given by (for a one-dimensional 
chain, similar expressions hold for higher dimensions) 



A = 2(u R u R ) 



h ■f—y COth^X; 



Nmui 



•?7 



B = 2(p R p R ) = — — 2_^xi COth^X;, 



(6) 



which arc independent of R because of translational sym- 
metry, and 



In this work we are interested in the states of thermal 
equilibrium, pp = j?e~P H , at inverse temperature f3 = 
1/ksT. For the harmonic lattice these can be written as 
a tensor product over all phononic modes, 



pp = 



(1 



■j—(3hijji 



(3) 



where hi are the number operators for each phononic 
mode. Since the Hamiltonian is quadratic the thermal 
states are Gaussian states (see for example [f| @ and ref- 
erences therein) which are specified uniquely by their first 



E = (urUp) + c.c. 



F = (prPp) + cc. 



h \ ^ / 2-kI \ coth^x/ 

> cos r , 

V \N J xi ' 



Nmuj 



mhio 
N 



V- /2ttZ \ 

2_^cos I — r J xi coth^x/, 



(7) 



which only depend on the distance between the two sites 
measured in terms of the lattice constant a, r = — — - G 
{1,2,3,...,^-^}. Here the inverse temperature is re- 
placed by a unit-free quantity £ = and the phononic 
frequencies are re-scaled by the interaction strength of 
nearest neighbours xi = —. 
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C. Two-site entanglement condition 

For Gaussian continuous variable states of two modes 
there exists a necessary and sufficient criterion that de- 
cides whether the modes are entangled. (The criterion 
is necessary and sufficient for Gaussian states of 1 x N 
modes, Q-) It is a variation of the positive partial trans- 
pose criterion (PPT-criterion) for discrete systems [1, 
For separability of two sites R and P in the lattice the 
criterion requires the positivity of a matrix inequality, 



Prp separable 4=> Y T Rp {p) + ihQ^a > 0, 



(8) 



R,P 



where a = 



1 

-1 



is the symplectic form. Tjfp is the 



covariance matrix with time-reversed mode P, 



1 RP ~ 



1 

-1 



RP 



1 

-1 



(9) 



i.e. for Tpp in Eq. ([5]) wc simply replace F by — F to 
obtain Tp^p. Using the symmetry of T^J, criterion © 
reduces to requiring the positivity of two temperature 
dependent functions, Si i2 (T, r), 



0<5i,2(T,r) = ^((u r ±up) 2 )((prtPp) 2 ) -1, 
n 



1 



(A±E(r))(BT F(r)) -1. (10) 



Here r = \R — P\/a is the separation of the two sites R 
and P in the one-dimensional chain. If positivity fails 
the two sites are entangled and the magnitude of their 
entanglement can be given by the logarithmic negativity 
E M , defined in [13, El- E N can be calculated directly 
from the functions Si 2 as 



Em = ^2 max{0, - In y/Sj + 1}. 



(11) 



Analytical and numerical results based on criterion ((Sj) 
for two sites are discussed in Section [TTT] and Section IIV1 
respectively. The criterion is also sufficient to detect en- 
tanglement between any two blocks of sites as shown in 
[l[. Another approach to reveal other kinds of entan- 
glement, i.e. exclude full separability, is the method of 
entanglement witnesses, introduced in 

[HE!. Entangle- 
ment witnesses are observables that can take a particular 
range of expectation values only for entangled state and 
therefore 'witness' entanglement. They have been used 
to investigate spins systems [H, [13, EH and also [l6j |. 
as well as Bosonic gases [l7| ■ In Section I VII we use the 
same strategy for the harmonic lattice to deduce a tem- 
perature below which entanglement must be present. A 
third approach to investigate entanglement in a harmonic 
lattice uses the full separability condition for the covari- 
ance matrix, T > (J^ Fj, for all sites j Q- The mathe- 
matical derivation following this direction is given in Q . 



Fig. |4] in the present paper unites the results on near- 
est neighbour entanglement and witnessed entanglement 
with these previous works and shows how the different 
strategies complete each other. 



III. TWO-SITE ENTANGLEMENT (ID) 

A. Entanglement in the ground state 

Any pair of nearest neighbours in the harmonic chain is 
coupled and we expect it to be entangled in the ground 
state. Before going in the discussion of the numerical 
plots in Section HVl let us here derive the negativity Eq. 
(|11[) for zero temperature and in the thermodynamical 
limit. 

Let T -> 0, i.e. /3 — > oo, then Si, 2 (r) = 
W f (1 ± cos ^ • r) (1 T cos ^ ■ r) - 1. Further 
substituting ^ = y and ^ = z and taking the contin- 
uum limit N — > oo to replace the sums by integrals one 
obtains 



Si, 2 W = 



1 ± cos 2zr 



(12) 



sin 



dyWsin y 



2^ 



(1 =F cos2yr) 1—1. 



Simplifying this expression further wc let 
S — * and by using the symmetry 
of the integrands one obtains Si ;2 (^) = 

dz i TS?ff ) (io f d y I sin yl (! T cos2yr)) 1, 
which can be easily evaluated for arbitrarily distant 
neighbours. For nearest neighbours we find 



(io 2 



Sx(r 
S 2 (r 



+oo, 
16 

37T 2 



- 1 ps -0.46. 



(13) 
(14) 



Finding S 2 (r = 1) negative implies that nearest neigh- 
bours are entangled and this can be quantified in terms 
of the negativity, Eq. (fTTj) . 



E N {r = 1) ps 0.31. 



(15) 



The behaviour of entanglement between any other pair 
of two sites in the chain is not obvious. Even though 
non-nearest neighbours do not interact directly, common 
nearest neighbours could mediate entanglement. Wc find 
that this is not the case when the chain size grows arbi- 
trarily, as both Si and S 2 remain positive for r = 2,3, .... 
However, for finite chain sizes, the above calculation is 
only an approximation and the numerical results show 
small entanglement between next nearest neighbours for 
very small chains at low temperatures, see Section ITVl for 
discussion. To conclude, at zero temperature and in the 
thermodynamic limit of large N two-site entanglement 
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only exists between nearest neighbours and not between 
any further separated sites. This result agrees with [is! ], 
where the authors showed that in gapped harmonic lat- 
tice systems the correlation functions decay exponentially 
in the ground state and thermal states. 



B. Critical temperature 

When the temperature rises, thermal mixing becomes 
stronger and correlations between the nearest neighbours 
decrease. We expect intuitively that any entanglement 
should vanish at large enough temperatures. The tran- 
sition into the classical regime is marked by a critical 
temperature which we here want to approximate. The 
exact temperature for various parameter ranges can be 
found numerically and is discussed in Sections HVl and Ivl 

Again we take the continuum limit N — > oo and 
replace the sums by integrals in Eq. (|10p . Assuming 
that the trapping potential vanishes, S — » 0, we ob- 

tain S 1)2 (r) = 4, (/* dz "> (1 ± cos2zr)) x 



J Q 2 dy coth (2£| sinyQ | siny| (1 qp cos2yr)J — 1. The in- 
tegrals have no closed expression and we need to further 
approximate the coth- function to evaluate Si^- The cou- 
pling parameter £ = where lo is the finite interac- 
tion strength between neighbours, is temperature depen- 
dent. Assuming high temperatures, T — * oo =>• (3 — > 0, 
one can expand to first order in £: coth (2£| sinzQ w 
1 Hil£££i, Thus for nearest neighbours one ob- 



2£| sinz| 

tains 



Si(r 
S 2 (r 



+oo, 
1 



< when £ > 1. 



(16) 
(17) 



Here we have dropped quadratic and higher powers of 
£ for consistency. Indeed 52 can become negative and 
the point where S^r = 1) crosses through zero gives the 
critical temperature for nearest neighbour entanglement, 
from £ nn . = 1 follows 



2kZ 



(18) 



We remark that this temperature is not as large as the ini- 
tial approximation assumed and one may wonder about 
its validity. Higher expansions could be considered to 
evaluate the critical value more accurately. We will not 
follow this direction here, instead we will present numer- 
ical data confirming the derived critical value in the next 
section. 

For the next-nearest neighbours the factors change 
slightly and the expressions become S\(r = 2) — > oo 
and S2(r = 2) = ^. Both expressions are never nega- 
tive irrespective of the temperature and no entanglement 
is expected between next-nearest neighbours in the con- 
tinuum limit. The numerical data in following section 
independently confirm this result. 




FIG. 1: The figure shows the negativity Em of a one- 
dimensional chain in thermal equilibrium over the temper- 
ature T (expressed in multiples of the coupling to). The trap- 
ping potential is small compared to the coupling, 8 — 10 _4 w. 
The solid curves show the nearest neighbour entanglement 
for an increasing number of sites: TV = 3 (uppermost), ...,51 
which converge to an asymptotic curve for N — > oo (low- 
est of the solid curves). The dashed curves show the next- 
nearest neighbour entanglement for an increasing number of 
sites, TV = 3 (uppermost), 15, which vanishes for further 
increasing N. 



IV. ENTANGLEMENT FOR FINITE 
TEMPERATURE AND SIZE (ID) 

A. Data 

Numerical data of nearest neighbour (solid) and next- 
nearest neighbour (dashed) entanglement in terms of 
the logarithmic negativity E_\f, Eq. lfTTj) arc displayed in 
Fig. [1] over the temperature and for varying size of a one- 
dimensional chain. We observe that in general the entan- 
glement is highest for low temperatures and decreases 
sharply with increasing T until it vanishes abruptly at 
some critical temperature. For nearest neighbours the 
critical temperature decreases for increasing TV and con- 
verges to T nn . ?s in the thermodynamical limit, 
which agrees with the analytical result in Eq. (fl8|) . Hence, 
the stronger the interaction u between nearest neigh- 
bours, the more robust their entanglement with respect 
to thermal noise. The maximal amount of entanglement 
at T = decreases for growing iV and converges to the 
threshold value of E^f ~ 0.3 independent of u, in good 
agreement with the analytical result Eq. P^|) . 

Next-nearest neighbour entanglement (dashed curves) 
is in general much smaller than the nearest neighbour en- 
tanglement. Nonetheless, the numerical data show that 
a small amount of entanglement exists between next- 
nearest neighbours when the chain is very small. How- 
ever, the entanglement as well as the critical temperature 
are considerably lower than for nearest neighbours and 
decrease rapidly to zero for increasing N. 
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B. Interpretation 

We can understand the behaviour of entanglement 
with changing temperature in the following way. With- 
out quantum zero-point fluctuations, the ions in the 
chain would be completely fixed and entanglement would 
stretch over all sites at zero temperature. However, due 
to the fluctuation of the phononic modes of the chain (# 
of phononic modes is N with N — ► oo) correlations are 
mixed with anti-correlations and the entanglement in the 
chain reduces considerably. For nearest neighbours the 
mixing of vibrations is biased. This is a result of the di- 
rect coupling between them leading to an average devia- 
tion of the ions from equilibrium position of less than the 
harmonic oscillator length of the interaction. Indeed for 
small temperatures entanglement exists between nearest 
neighbours, see the solid curves in Fig. Q] 

Next-nearest neighbours and more distant sites do not 
have a bias of the phononic modes in the thermodynamic 
limit and entanglement between them is averaged out 
even at T = 0. However, when the chain size is small, 
finite size effects occur. So few phononic modes are avail- 
able that the different swinging motions of the ions are 
not completely balanced and some correlations survive. 
This is why we find next-nearest neighbour entanglement 
for very small lattices, see the dashed curves in Fig. [T] It 
is also the reason why the entanglement between nearest 
neighbours is in general slightly higher for smaller lat- 
tices, as can be seen for the solid set of curves where the 
uppermost curve displays the behaviour of a tiny chain 
with N = 3, and the lower curves apply for increasing N. 

When the temperature is raised, more phonons are ex- 
cited and the mixing of different vibrations is enhanced. 
This leads to a rapid decay of the entanglement and when 
the temperature reaches the magnitude of the phononic 
modes that drive nearest neighbours against each other 
the entanglement vanishes completely. This happens 
when the two degrees of freedom, the position and the 
momentum of each of the iV ions, with classical thermal 
energy per degree of freedom exceeds the energy 
of the quantum fluctuations of each of the ./V springs 
that connect nearest neighbours, %^ (for 8 f=a 0). We 
thus expect nearest neighbours to behave classically when 
2 , k^r_ ^ nw j g fulfil]^ This heuristic argument gives 

the critical temperature for entanglement, T n . n . = jj^-, 
in agreement with the numerical cutting point in Fig. [l] 
It also indicates that the energy is a good witness of en- 
tanglement, which wc will elaborate in Section [VII 

V. DEPENDENCY ON FREQUENCY 
SPECTRUM AND DIMENSION OF SPACE 

A. Data 

The impact of the phononic frequency spectrum on 
the entanglement properties are studied quantitatively 
by varying the ratio between trapping potential and 



0.3 ~ 



0.2 



0.1 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2k B T 

FIG. 2: Entanglement, Em, of nearest neighbours in a one- 
dimensional chain with N — 49 sites, depicted over tem- 
perature, T. The different curves represent increasing ra- 
tios of the trapping potential versus the coupling strength, 
with the uppermost curve reaching Em m 0.3 at T = for 
8/lu — 10~ 4 , and lower curves for the parameters 5/oj = 
10" 7/2 , 10~ 1/2 , 1 (dashed), Vw . 



nearest neighbour coupling. The numerical results arc 
displayed in Fig. [2j As before, for small 8/u> the en- 
tanglement starts at Em ~ 0.3 for T = and de- 
creases for increasing temperature until it vanishes at 
7n. n .(£/w << 1) « However, when the trapping 

potential is raised and reaches the magnitude of the in- 
teraction the initial entanglement starts off much lower, 
Em ~ 0.17 for 8 = u>. Notable is that this entanglement 
is more robust and survives increasingly high temper- 
ature mixing. For example, when 8 = u>, the critical 
temperature has increased to T n n .(i5 = us) « 1.25^-. 



B. Interpretation 

When the trapping potential 8 is raised close to and 
above to two effects emerge. Firstly, the on-site potential 
fixes the ions tightly to their sites and only a small frac- 
tion of the available energy can be put into the interac- 
tion. Even though they arc coupled, the ions cannot fol- 
low the motion of their neighbours and entanglement be- 
tween nearest neighbours starts off at substantially lower 
values for T = compared to the case of small 8, see 
Fig. [21 The second effect is that higher phononic modes 
'cost' a higher amount of energy. Thus for growing 8 the 
bias of phononic modes cannot be compensated even at 
higher temperatures and nearest neighbours stay entan- 
gled up to increasing critical temperatures, as shown in 

Fig. m 

For higher dimensions, that is, when the chain becomes 
a two- or three-dimensional square lattice, the behaviour 
of two-site entanglement is qualitatively similar to the 
one-dimensional case, as shown in Fig. [3] Comparing 
this figure with the graphs for varying 8 in Fig. [2] one 
can see that raising the dimension leads to similar effects 




FIG. 3: Entanglement of nearest neighbours in ID (solid line, 
N = 31), 2D (dashed, N = 31 2 ) and 3D (points, N = 31 3 ) 
is shown over the temperature. The trapping potential is 
6/cj = 1CT 4 in all cases. For a discussion please refer to the 
main text. The dimension- independent temperature for any 
witnessed entanglement, Tew, derived in Section IVII is also 
indicated and discussed in Section TV HI 



as raising the trapping potential. The higher the dimen- 
sion the lower is the ground state entanglement, but the 
higher rises the critical temperature. 

The reason behind this behaviour is the growing fix- 
ation of the ions to their sites for increasing dimension. 
Namely, in the thermodynamic limit (and for S — ► 0) the 
average displacement of the ions from their equilibrium 

- 2 

position, (us ), diverges linearly in ID and logarithmi- 
cally in 2D, while converging in 3D. No long-range order 
can be established in ID and 2D as the crystal simply 
'melts' (for a discussion see for example [|J and references 
therein) . Only in three dimensions do the fluctuations of 
the ionic positions become smaller than the lattice spac- 
ing, long-range order can occur and a solid with a three- 
dimensional crystalline structure can form. To obtain 
the same effect one could impose a strong trapping po- 
tential to hold the ions tightly in place and thereby force 
crystallisation in low dimensions. Thus, increasing the 
trapping potential is effectively raising the dimension of 
the lattice structure and vice versa. 



VI. ENERGY AS AN ENTANGLEMENT 
WITNESS IN d DIMENSIONS 

Two-site entanglement is a very restricted kind of en- 
tanglement. Other kinds of entanglement exist, for exam- 
ple genuine three-partite entanglement. As indicated by 
the previous numerical results the balance between the 
available thermal energy compared to the phononic spec- 
trum determines the existence of entanglement implying 
that the energy may be a useful criterion as an entan- 
glement witness. This idea has been used in a number 
of works before, for example in [l3|, 0, HI, E3]- Adopt- 
ing this strategy for the harmonic lattice in d dimensions 



we identify the temperature range where some kind of 
entanglement must necessarily be present. 



A. Derivation of energy bound for separated sites. 

The most general, fully separable (into N sites) state 
is of the form 



Pa 



E 



(19) 



where pj > arc probabilities with pj = 1 . The en- 
ergy (H) (c.f. Eq. ([I}) can be evaluated for all separable 
states as 



{H) sep = J^Pj tr 



R I p\ md 2 u 2 R 2mJ 2 u 2 R 



2m 



-mw 2 Pi tr [Pf (««)] tr [Pf ("«')J ■ ( 2 °) 

The second term can w.l.o.g. be set zero since all first 
order expectation values of any operator can be shifted 
locally to zero without affecting the non-local entangle- 
ment properties. Thus, the energy of the whole lattice 
in any separable configuration, p sep , becomes the sum of 
the individual energies for all sites, 



(H)sex> = 



5> 

R 



P R H R 



(21) 



where p R = ■ pj p R is an arbitrary state for site R 
and^:=H + ^ + . 



^ the effective Hamilton 



2m 1 2 1 2 
operator for one site R alone, where the ur and pn are 

in general d-dimensional operators. Thus the description 
of the lattice has become a mean-field model where each 
site interacts with an effective, identical background pro- 
duced by the interaction with all other sites. Reversely, 
a description using only effective one-site Hamiltonians 
leads to a classical model where entanglement has been 
averaged out and states are always separable. 

H R can be rewritten as a harmonic oscillator, H R = 
TiSl (n R + |) , where the frequency f2 = \/2uj 2 + S 2 is 
identical for all sites R. By the presence of quantum 
zero-point fluctuations we always have (Hr) > | hQ for 
any state p R . The possible energy for a separable config- 
uration of the whole lattice is thus bounded from below 
by 



(ff) 8 ep-sites > 2 m = 2 NhQ ' 
R. 



(22) 



Any state having an energy below this bound is neces- 
sarily entangled with respect to the sites. 



B. Entanglement for different modes. 

It may seem puzzling that one could apply the same 
strategy for other set of modes, for instance, for the 
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phononic modes themselves and find a non-zero energy 
bound. Indeed, for the normal modes one would find 

(-ff)scp-phonons > ^ (^) 
F 

with as given in Eq. |(5J) for <i dimensions. This in- 
equality is always fulfilled as the right hand side gives 
just the ground state energy. The criterion is thus use- 
less here as it can never be violated and is hence unable 
to detect any entanglement. 

By symmetry we expect that the sites in space are 
'most' entangled, i.e. give the highest energy bound, 
while the phononic modes remain separable in any ther- 
mal state. The situation for 'intermediate modes' is 
between these extremes. To show this we use again 
the energy witness argument (assuming one dimension 
for simplicity). For all orthogonal transformations O 
that transform the set of phononic momentum opera- 
tors, Vi, into some new momentum operators Pk, such 
that Pk = 53/ Oki Vi and similarly for the position oper- 
ators, Uk = 53; ®ki U\, so that the commutation relations 
remain [Uk,P m ] = ifiSkm, the Hamiltonian becomes 

~ ^(Pf m^U?\ yr^fPl mWl k Ul\ 
H = ^{^ + ^)=^{^ + ^t^) 

+ V- rnWl m U k Um (24) 

k^m 

where the interactions are given by the elements of the 
matrix W 2 m = 53; ®ki Of m . With the same argument 
as before we obtain a minimal energy threshold for all 
separable states with respect to the new set of modes, 

<ff) 8 ep-fc- m odes > \ £ h^W^. (25) 
fc 

In particular, the diagonal frequencies are Wf fc = 
53; with the doubly-stochastic matrix O 2 , 

5], 2 kl = 1 = 53 fc 0\ v With the Birkhoff-von Neu- 
mann theorem O 2 can be written as a convex combi- 
nation of permutation matrices of the same order and 
because the square-root is concave we only need to con- 
sider the two extremal cases which bound 53fc V^Ife 
from below and above. The first is a single permu- 
tation n with 0\\ = 5^n(k) resulting in a reordering 
W 2 k = w n(fc)' an< ^ the second is the complete mixture 
of all possible permutations such that 0\ x — jj and 
Wkk = W Si w f = These extremes are indeed iden- 
tical to the phononic modes and the sites, respectively. 
The energy bound for the general set of modes is thus 
bounded from below and above, 

being bigger than the zero-point energy, while remaining 
smaller than the energy bound for individual sites. Site 
entanglement thus possesses the highest energy bound. 



C. Critical temperature for separated sites. 

The energy of the d-dimensional harmonic lattice in a 
thermal state can be related to the temperature. The in- 
ternal energy, U, can be explicitly calculated and for high 
temperatures one finds the Boltzmann's equi-partition 
law, where each degree of freedom of the system obtains 
a thermal energy contribution of ^aX 5 

U(T) d Nk B T for T -> oo. (27) 

(In this context the relation is identical to the Dulong- 
Petit law of constant heat capacity at high tempera- 
tures.) Setting the internal energy equal to the energy 
bound for the sites, Eq. (|22|). one finds that any fully sep- 
arable configuration has a temperature above a certain 
critical value Tew 

fiQ Hlu I S 2 , . 

t -^ = 2^V 2 + 3^ Tew - (28) 

Conversely, any configuration with a lower temperature 
must necessarily contain entanglement between the sites 
of the lattice. Remarkably, the critical value is inde- 
pendent of the dimension of the system. Furthermore, 
the witness temperature for S = 0, Tew = 4g ~ 

1.4142 jjt^ is close to the nearest neighbour critical tem- 
perature T n n . w in ID. In higher dimensions the 
nearest neighbour critical temperature increases and ap- 
proaches the witness-temperature, as shown in Fig. [3l 

The energy relation Eq. (|27[) and hence the witness 
temperature Eq. (|28|) is initially only valid in the asymp- 
totic limit of high temperatures. For small temperatures 
the behaviour of the internal energy U (T) and its deriva- 
tive, the heat capacity CV(T), has been subject of intense 
discussion at the beginning of the 20th century (which 
is covered in many textbooks, see e.g. [H|). Introduc- 
ing a heuristic cut-off frequency lod, Debye found that 
in three dimensions the phononic oscillations lead to a 
heat capacity proportional to T 3 , implying Ud oc T for 
temperatures below Tn = tt 2 -. With the same method 

the internal energy scales with Ud oc T d+1 in lower di- 
mensions. Discussing the Debye temperature in detail is 
beyond the scope of this paper, however, direct conclu- 
sions can be drawn just from the properties of the energy 
function. Because of the monotonicity of the internal en- 
ergy and its convexity at small temperatures, Ud(T) is 
below the value of the high-temperature limit, U(T) in 
Eq. I[27p. in the regime T < Tq. The witness tempera- 
ture derived with the high temperature approximation, 
Tew, is thus a lower bound to the critical temperature 
that one would obtain when using the low temperature 
formula. The entangled region detected with the high 
temperature approximation is thus a subset of the region 
which one could detect when using the small temperature 
approximation. This implies that Tew is a valid bound 
for witnessed entanglement also at small temperatures, 
being most tight in one dimension and lesser so in higher 
dimension due to the stronger bend of the T d+1 -law. 
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FIG. 4: The graphic keeps the exact scale relations of the 
different phases of entanglement and separability in the har- 
monic chain (iV = 49), for varying trapping potential 8 and 
over the temperature T. Nearest neighbour entanglement ex- 
ists in the dark shaded region up to the temperature thresh- 
old T n . n .. Entanglement is witnessed by the energy in the 
hatched area, with threshold temperature Tew- Additionally 
the temperature threshold, Tbiocks (black point), for the exis- 
tence of entanglement between two neighbouring blocks of size 
N/2 from [lj (for a single value of 8) is displayed. The chain 
passes into the light shaded region of complete separability at 
the temperature T cr i t from [j| • For 8 = the threshold values 
are displayed below the figure, for their discussion please refer 
to the text. 



VII. ENTANGLEMENT PHASE DIAGRAM 

Fig. [4] is a comprehensive and exact diagram com- 
paring the two derived temperature thresholds for near- 
est neighbour entanglement and witnessed entanglement 
with previous results from and 0. As found in Sec- 
tion IIVI two-site entanglement resides between nearest 
neighbours (dark shaded) for temperatures below T n . n . w 
0.5 TiLo/kB (for N — > oo, 6 = 0) while non-nearest neigh- 
bour entanglement vanishes for large N. The hatched 
set are states whose non-separability is detected by the 
energy witness. It fully includes all nearest neighbour 
entanglement and extends further up to the temperature 
Tew ~ fitd/kB, cf. Eq. (|28[). where no two-site entan- 
glement exists. A different phase of how entanglement is 
distributed is thus reached. 

For example, entanglement can exist between two 
neighbouring blocks of N/2 sites. Results in Q for <5 = 
u> I V20 (this appears to be the critical value according to 
Figure 12 in [l| after retaining all units) confirm this 
and their critical temperature for block-entanglement, 
^blocks ~ 0.63?kj/fcB, is indicated in Fig. [4] by the dark 
point. When the temperature is raised further the chain 
passes over into complete separability (light shaded), i.e. 
no entanglement exists for any possible split. The tem- 



perature at which this transition happens is given as 
T crit = hu/1.2k B ~ 0.83hcu/k B (for 5 = ) in [|. 
For sufficiently small 6, the witness bound Tew is in 
good agreement with the exact threshold for full separa- 
bility, T cr it [H . This indicates that the temperature range 
where the entanglement survives thermal mixing is iden- 
tical to the regime where the quantum zero-point effects 
are important. The use of the energy as an entangle- 
ment witness is thus justified and in general an alterna- 
tive, easier way to evaluate regions where entanglement 
is present, cf. [ID, Q EE E3 ■ For higher dimension we 
find that the critical temperature for nearest neighbour 
entanglement increases towards the dimension indepen- 
dent witness temperature (see Fig. [3]) thus leaving less 
temperature range where non-two-site entanglement is 
witnessed. 

When 8 is large compared to uj, the frequency spectrum 
shrinks to essentially one dominating frequency, u>i w S. 
The mixing of frequencies occurring at increasing temper- 
atures becomes irrelevant and entanglement persist up to 
higher temperatures, the higher 6 becomes. The asymp- 
totic behaviour of the witness-temperature becomes pro- 
portional to S as kB j[™ oc ^j. Yet the full-separability 
temperature scales with ^ as fc ^J° ~ 2 which i m ~ 

plies a growing gap between the witnessed and the true 
entangled region for growing 8. 

Interestingly in this model any entanglement must be 
non-PPT [20J, in other words, there always exists a bi- 
partite split of the chain into two parties such, that par- 
tial transposition on one of the parties results in a non- 
positive density matrix. Any entanglement can thus po- 
tentially be identified by this property and the value 
of the negative eigenvalue can be used to measure the 
amount of entanglement. 

VIII. CONCLUSION AND OUTLOOK 

In condensed matter physics the two-point correlation 
function plj | between two sites is an important quantity 
used to determine the behaviour of the material. If cor- 
relations stretch over a wide range of sites, i.e. the cor- 
relation function becomes a constant or decreases poly- 
nomially for increasing distances, then one speaks of the 
establishment of an ordered phase. Quantum physics al- 
lows for a stronger kind of correlation, that is, entangle- 
ment. The entanglement criterion Eq. (| 1 0[) is reminiscent 
of the classical order criterion for the two-point correla- 
tion function. Instead of requiring that the fluctuations 
of the ionic positions, (u%), should be smaller than the 
lattice spacing, we here have a condition on the prod- 
uct of the relative positions and the relative momentum 
fluctuations. In other words, instead of requiring order 
in space the entanglement criterion (|10[) requires order 
between two sites in phase space. 

We have found that entanglement exists between near- 
est neighbours in all dimensions, for all lattice sizes, and 
for all values of S, when the system is cooled below some 
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critical temperature, which depends on the former pa- 
rameters. By changing the temperature the entangle- 
ment can be modulated and even cut out. The critical 
temperature at which the entanglement vanishes is di- 
rectly proportional to the coupling constant; for the har- 
monic chain this is the coupling lu. This relation is in 
agreement with previous results in spin chains [l3| and 
Bose gases [l7l . l23j where also the interaction strength 
determines the critical temperature. 

Entanglement between next-nearest neighbours exists 
in very small lattices, but is in general much weaker in 
magnitude and also only observable at much lower tem- 
peratures than for the nearest neighbour case. A possible 
application for trapped ions in lattices with harmonic in- 
teraction, would be to tune the value of S and T such, 
that only nearest neighbour ions remain entangled. By 
slowly inserting an additional ion in the lattice one can 
destroy entanglement to the former nearest neighbours 
and the new nearest neighbours become entangled. Vary- 
ing the trapping potential 5 or the interaction strength ui, 
one can further change the amount and range of entangle- 
ment in the lattice. Such techniques and the possibility 
to apply additional external potentials that modify the 
interaction and hence affect the structure of the entan- 
glement may allow to store and manipulate information 
on the ions. For instance, applying an extra harmonic 
potential to all ions in a small periodic chain results in a 
higher frequency of the lowest phononic mode, the centre 
of mass mode, that could influence the chain to 'swap' 
entanglement from nearest neighbours to opposite ions 
in the ring by simply raising the temperature. 



Even though in the thermodynamic limit only neigh- 
bouring pairs are entangled and more distant pairs are 
not, we can speak of an 'entangled phase'. This notion 
is justified because entanglement is much stronger than 
classical correlations and nearest neighbour entanglement 
is already sufficient to transport information across the 
whole lattice. Such correlations enable universal quan- 
tum computation in the measurement-based model, as 
initiated in (24j for spin systems with nearest neighbour 
interaction and its analogue for the continuous variable 
case proposed in [25[. The entangled phase survives, ap- 
proximately, as long as the system is in the quantum 
regime where the zero-point fluctuations dominate the 
physics of the lattice. Going beyond that range implies 
the validity of a mean-field theory where each site can 
be described as an independent oscillator whose effective 
frequency includes the contributions from the rest of the 
system. 
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